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Abstract 

By solving an infinite nonlinear system of g-difference equations one constructs a 
chain of q-difference operators. The eigenproblems for the chain are solved and some 
applications, including the one related to q-Hahn orthogonal polynomials, are discussed. 
It is shown that in the limit q ^ 1 the present method corresponds to the one developed 
by Infeld and Hull. 
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1 Introduction 

The discretization of the ordinary differential equations is an important and necessary step 
toward finding their numerical solutions. In place of the standard discretization based on the 
arithmetic progression, one can use a not less efficient g-discretization related to geometric 
progression. This alternative method leads to g-difference equations, which in the limit q ^ 1 
correspond to the original differential equations. The theory of g-difference equations and 
the related g-special functions theory have a long history (see e.g. (Ij). During the last two 
decades they have been reviewed because of the great success of the theory of quantum groups. 

The other crucial way of solving ordinary differential equations is based on the factor- 
ization method first used by Darboux [2]. Later the method was rediscovered many times, 
in particular by the founders of quantum mechanics, see [8] and [11], while studying the 
Schrodinger equation. We refer to [16] for an exhaustive presentation of the factorization 
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method. In the paper [2], which is now considered to be fundamental, Infeld and Hull sum- 
marised the quantum mechanical applications of the method. Fixing an infinite system of 
Riccati type equations they have constructed a chain of second order differential operators 
and proposed some method of solving corresponding eigenproblems. 

In this paper we construct the chain H2.71|) of second order g-difference operators by solving 
an infinite nonlinear g-difference system. This chain depends on a freely chosen function and 
a finite number of real parameters. In Section 2 we find a family of eigenvectors for the 
operators of H2.7H1 . In Section 3 it is shown that g-Hahn orthogonal polynomials, which are 
g-deformation of the classical orthogonal polynomials, form the family of solutions obtained 
by our method. Other examples of solutions obtained by the factorization of g-difference 
equations are presented in Section 4. Finally passing to the limit g ^ 1 in H3.36|l . (|3.37p we 
obtain some new families of solutions for second order differential equations. 



2 Factorized chain of the second order ^'-difference operators 

In this section we shall consider the sequence of the second order g-difference operators 

Uk = Zk{x)dgQ-^dg + Wk{x)d, + Vkix) , keNU{0} (2.1) 

acting in the Hilbert spaces Tik- By definition 7ik consist of the complex valued functions 
■0 : [a, 6]q — > C defined on the g-interval 

[o, b]g := {q^'a : n G N U {0}} U {q% : n E N U {0}} (2.2) 

and square-integrable, i.e. {tp\ip)k < +oo, with respect to the scalar products 



{'il^\f)k ■■= J ip{x)(p{x)0k{x)dqx . (2.3) 

[a,b]q 

Let us recall (see [34]) that by definition the g-derivative is 

(1 - q)x 

and the g-integral on the g-interval [a, b]g is given by 

/oo 

i>{x)dgx := ^(1 - g)g" (5V'(g"6) - aV'(g"a)) • (2.5) 

la,b\q 

If o = and b = oo then 

/ ^{x)dgx := lim / ^{x)dgx = V (1 - • (2.6) 

/ n—>oo I '—^ 

"=-°° 

In the case if a = — oo and b = 

f ^P{x)dgx := lirn^ I ij{x)dgx= (1 - g)g" (/(g") + /(-g")) . (2.7) 
In the limit q —> 1 the above definitions correspond to their counterparts in standard calculus. 
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The scalar products (|2.H|1 are defined by the weight functions Qk '■ [a, b]g — > M, which are 
related by the recursion relations 

Qk-i = VkQk (2.8) 

and 

Qk-i = Q{Bkgk) , (2.9) 
where ry^, Bk are real valued functions on [a, b]g and the operator Q is defined by the formula 

Qf{x) = (p{qx) . (2.10) 
For the sake of consistency we need to add the conditions 

Q{BkQk)=VkQk (2.11) 

on the functions r]k and Bk- Additionally we impose the boundary conditions 

Bkia)Qk{a) = Bkib)Qkib) = . (2.12) 

If we introduce the functions 

we can rewrite the formula H2.1H1 in the form of a g-Pearson equation jl,S| 

dg{BkQk) = AkQk. (2.14) 

In the limit g — > 1, the equation H2.14jl corresponds to the Pearson equation which is important 
for the theory of classical orthogonal polynomials 

We say that the operators admit a factorization if 

Uk = AlAk + ak (2.15) 

and 

Hfc = dfc^ jAfc+i A*+i + Ofc+i) , (2.16) 
where the annihilation operators : Tik 'Hk-i are of the form 

Ak = dg + fk (2.17) 

and fk are real valued functions on the set [a,6]g. The adjoint operators A^ : Tik-i Ti-k, 
called the creation operators, are given by 

= {dg + fkT = Bk {-dgQ-^ + fk) - Ak (1 + (1 - q)xfk) . (2.18) 

The derivation of the formula (|2.18|1 is given in Appendix A. It follows from l|2.15|l that the 
real valued functions Zk-, Wk and Vk are related to fk, Bk, Ak by the formulas: 

Zk = -BkQ-^ (1 + (1 - q)id fk) , (2.19) 
Wk = Bkfk - Afc (1 + (1 - q)id fk) - q-^BkQ-\fk) , (2.20) 
Vk = -Bkdg {Q-Hfk)) - Akfk (1 + (1 - q)id fk) + + . (2.21) 
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Necessary and sufficient conditions for the consistency of factorization formulas (|2.15|l and 
itTTTni are 

rjk+iix) = gk{x)rik{q^^x) , (2.22) 

ipk+i{x) = ^^^^'fkiq'^x) , (2.23) 
9k[x) 

Mx) - ^^akiqx) = (2.24) 

q^dk+iBk{qx) - gk{q^x)Bk{q^x) \ gk{qx) 



{l-qfq^x^ + dk+iak - ak+i J 

where we have introduced the additional notations 

^ . (2.25, 

ifkix) ■■= fk{x) + ,^ ^ . , (2.26) 
Okix) := ipl{x)7]k{x) . (2.27) 

The detailed derivation of these formulas is given in Appendix B and in the paper jT2l l5]. 

The relations (|2.22p . (|2.2,Sp and (|2.25p . (|2.27p allow us to express the functions B^, rjk, ^k 
and Ofc by the initial data i?o, and oq 

Bk{x) = gk^i{x)gk^2{x) ■ . . goix)Bo{x) , (2.28) 
Vkix) = gk-i{x)gk-2iq'^x) . . . goiq'^^^x)r]oiq'''x) , (2.29) 

Substituting H2.28M2?8T1l into condition (j2.24p we obtain the infinite sequence of the non- 
linear functional equations 

ao{x) - dk+i^^^^j-^ao{qx) = Gk+i{x) {dk+iau - afc+i+ (2.32) 
Gk{qx) 

^ q^dk+igk-ijq^^^x) . . . go{q^+^x)Bo{q^+^x) - gfc(g^+^x) . . . go(g^+^x)^o(g^+^x) 

(1 — q)"^ q'^'^'^^ x"^ 

where 



gk- iiq''x) . ■■go jqx) 

(dk-.-diY 

Goix) := 1 , (2.34) 



Gk{x) := ^"-'Y^ -'■■■^^y'^-) for A:GN, (2.33) 



for the functions oq, Bq and gfc for A; G N U {0}. 

One sees from (|2.28M^l?T1l that the sequence of functions gk, k S N, satisfying H2.32p 
defines the chain of q- difference operators (j2.ip if the first element Hq of the chain is given. 
So, the problem of construction of the factorized chain given by H2.15p and (|2.16p is equivalent 



to solving of the system of functional equations (|2.32l 
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Let us now present the limit behaviour of the formulas obtained above when the parameter 
q tends to 1. It is easy to see that the set [a, b]g becomes the interval [a, b] in the limit g — > 1 
and the scalar product turns to be 



{■ip\ip)k = / ■4'ix)Lp{x)gkix)dx , (2.35) 



where the weight function Qkix) satisfies the Pearson equation 

-j^iQkBk) = QkAk , (2.36) 

with the boundary conditions (|2.12p . For g — > 1 the operator Q goes to the identity operator 
and d„ > 4^. In the limiting case the annihilation and creation operators are of the form 

Afe = ^ + fk, (2.37) 

Al = Bk(^-^ + fk^-Ak (2.38) 
and the operators are given by 

^'^ " ~ ^^1^ + " ■^^^^'^ " ^^^^ + ■ ^^-^^^ 

The g-difference equation (|2.H1 tends to the differential equation 

(P d \ 

^'^^^^d^ + ^'^^^^d^ + ^^^""^J "^^^""^ " ^kMx) , (2.40) 

where the coefficients are given by 

Zkix) = -Bkix) , (2.41) 
Wk{x) = -Ak{x) , (2.42) 
Vk{x) = {fiix) - /^(x)) Bk{x) - fk{x)Ak{x) + ak . (2.43) 

The recurrence transformations H2.22MT^r?|l for g — > 1 tend to 

Bk+i = dk+iBk , (2.44) 

Ak+i = dk+i {Ak - -^B^ ■ (2-45) 

The sequence of g-difference equations H2.24|l tends to the sequence of non-linear differential 
equations 

Bk{fl+i - fk + f'k+i + fk) - Mfk+i - fk) + 2B'jk+i - K + B'l = ak-'^, (2.46) 

k G NU{0}. The equation (|2.46p for Bk{x) = 1 and Ak{x) = was considered in many papers 
(see [3 El El E» ilSj), but nevertheless for these differential-difference equations there is 
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no complete theory. One of the methods for solving of (|2.4()|l is to look for the solutions of 
II2.46I1 in the form of infinite series 



A = i;/*(x)fc* (2.47) 

and obtain in this way the conditions on the function fi{x). The case of solutions given by 
the finite series were consider by Infeld and Hull [Jj. The classification of all factorisable 
one-dimensional problems is still an open question. 

Now, we come back to the general case. Regarding the extreme nonlinearity of the system 



H2.32|l . the possibility to solve it is rather out of the question. Therefore, we shall restrict 
ourselves to the subcase 

gk{x) := dk+iq^ for 7 e M (2.48) 
and consider the system (|2.H2p . which is reduced now to 

ao{x) - q'^aoiqx) = — (4+iafc - Ofc+J + (2.49) 

«fc+i . . . «i 

I «2(fc+i)7ofc+i <l^~^Bo{x) - Bo{qx) 
(1 — q)'^qx'^ 

as the infinite system of equations on the initial functions Bq and oq. Eliminating oq from 
H2.49|l we obtain 

(1 - qfq^^'^d^^x'^ — {dk+ittk - at+i) - ditto + ai] = (2.50) 

V"fc+i..."i / 

= q'-'Boiqx) - Bo{q\) - q^^^^''^^ [q^'^ Bo{q'^+' x) - Bo{q^+^x)) , A; G N . 
Now, we shall look for the solution of H2.5()|l in the form 

So(x) =x'^^6„x" , (2.51) 

where 5 € (0, 1). Substituting H2.5H1 into H2.5()|l and comparing the coefficients in front of x"^ 
we obtain the expressions for the G M 

afe+i = 4+1 . . . diq-^^ (-"0^ + ^^^^7^ " «^2[7fc][7(fc + 1)) , A; G N (2.52) 
and the function Bq 

Bo{x) = b2X^ + + box^-^^ , (2.53) 
where 62,^*1,^0 G ^- From ^Hu^ and IfT^ we have: 
(i) if 7 / 0, then 

f X q'^^'^h q"'{ditto -ai) ^^"'^^0 -27 .0 ka\ 

where /i G M; 
(ii) if 7 = 0, then 

ao{x) = h and dioo = cli , (2.55) 

where /i G M. 
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Finally, substituting H2.48|l to H2.28MT^TT|l we find the following transformation formulas 

Bk{x)=q^''dk...diBoix) , (2.56) 

mix) = q^'^'dk . . . diTjoiq-'^x) , (2.57) 

Mx) = q-^^'MQ-'x) , (2.58) 

afc(x) = q-"'''dk . . . diao{q-''x) , (2.59) 

where Bq, are given by (j2.5H|l and (|2.54MTKK|) respectively. The functions % and ^Poix) 
are related to and ao by 

Vo{x) = b2x'^ + + box^-'^^ - (1 - q)xAo{x) , (2.60) 

At the moment, given the functions Bq, ao, we can use (|2.56MT59)l . H2.13Mm|) . (|2.26p and 
H2.27|l in order to express the functions A^, and Qk'. 



Ak{x) =q'<^dk ...di (q'^Ao{q-^x) + [-2k]h2X+ (2.62) 
+ [A:(7 - 2)]hix^-^ + [2fc(7 - l)]hox^'^^) , 

fk{x) =q-^^Uq-^x) - ~ , (2.63) 

(1 - q)x 

/ N g ^ go(g a;) 

^^■^"^^ ... 4 nn=0 (&2g-2"x2 + 6ig«(7-2)3;2-7 + boqMl-1) x^-^"t) 

by Ao, /o and qq. From conditions (|2.13Mm|l . H2.26|l and (|2.27p we see that the functions 
^0, /o, ^?o are related by 

= + + 60X^-2, _(i_,)^„(,)^>o(^x) , (2.65) 

1 \ ^ Cto(3;) 
^l_g)x) " 62^2 + 6ix2-7 + 6ox2-27 _ (1 _ g)xAo(x) ' ^^'^^^ 

So, further we shall assume that the function ^^|gy is continuous in 0. Under this assumption 
we obtain from l[T65ll and (l2^ 



■^"^^^ V 62X2 + 61X2-7 + 6ox2-27- (l-g)xAo(x) (1 - q)x ' ^^'^^^ 



oo 



This means that one finds the explicit formulas for the annihilation and creation operators 



=4 . . . di (-^^'=(62x2 + 6ix2-^ + 6ox2-2->) (^S.Q-^ + (^-^) + 

+ Y''"o(g"''x)r/o(Q'"''x) 



(2.70) 
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and from this the explicit expression for H^, 

Hk = dk ■ . . di -(1 - q)q~^x^{b2 + bix~^ + b^x' 



dg+ (2.71) 



+ (-q-'x\b2 + bix"' + 5ox-2^)W 2{!i-t+'^)x) + ^"o{q-''x)r^o{q-'x)\ dg+ 



+ 



b2 + bix-"/ + bpx-'^"' 
(l-qf 



q- {I- q)x\ 



^ao(g-(fc+i)x) 
rio{q~^''~^^^x) 



+ 



+q ^^Mq~ 



—^\/vo{q ^x)ao{q ''x) + 



(1 — q)x 



-Mk-l) [ [7(^-1)] Ol M , , r ,1 

-9 ^(^ao -^ + <762[7(A:-l)][7fc] 

which depend only on a function Aq. 

The chains of operators A^., and appearing in H2.69p . (|2.7()p and H2.7H1 in the limit 
g ^ 1 are given by 

Afc = — + /o(x) H , 



A^ = 4 . . . c^i {bo{x) {-^ + Mx) + ^^^^^) - Mx) + k^Boix] 

/ d 
Hfc = 4 . . . di f -So(x)^ - (^o(a;) - kB'Q{x))—+ 



where 



/o(x) + (^Aoix) - kB'oix)) - ao{k -l) + ^k- b2l^k{k - 1) 

XI d\ 



Bkix) = dk ■ ..diBoix) 



(2.72) 
(2.73) 
(2.74) 



(2.75) 



Akix) =dk... 4(^0(2;) - k—Bo{x)) 



fk{x) = fo{x) + k{^ - 1)- 



Qk{x) 



Qq{x) 



dkdl_^...d\ B^{x) 
and the functions i?o, /o and qq have the form 



(2.76) 
(2.77) 
(2.78) 

(2.79) 



Mx) 



-b2{-f+l)x+^:^^^i^^x-bihx^-'' -bo{l-y)x^-^-' +Aoix) 

^bixHbix'^-'y+box'^-'-i'y) dla 7 7^ 

dla 7 = 



al _|_ Ao{x) 
' 2x^ 2(fe2+bi+bo)3;^ 



(2.80) 
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1 Mllrlt 

Qo{x) = ■ (2.81) 

Summing up we see that the construction presented above gives us the nontrivial chain of 
Hamiltonians H2.7H1 parameterised by the freely chosen function Aq and the real parameters 
60, bi, 62, h, dk and 7, G N U {0}. 



3 Eigenvalue problem for the chain of operators 

We shall be interested in solving the eigenvalue problems 

HfcV'fc = AfeVfc for keNU {0} . (3.1) 

If the operators H^, admit the factorization given by H2.15p and H2.16|l then the eigenvalue 
equation (|3.H) is equivalent to the two equations 

AlAkipk = (Afc - ak)tpk , (3.2) 
Afc+iA^+iVfe = (4+iAa; - afe+OV'fc • (3.3) 

From and (pO|l one gets 

iikAl_^i'4>k = dk+iakAl_^i'4>k (3.4) 



if 

or equivalently, if 

Let us remark here that 



Hkipk = aki^k (3.5) 
AfcVfc = . (3.6) 

Afc+iV'fc = (dfc+iAfc - Ofc+i) V'fc (3.7) 



and thus A\j^^il^k £ "^fc+i- The formulas (|3.7p show also that the application of A^+i to 
A^_l_^Vfc turns it back to the eigenvector of Hfc proportional to the eigenvector ipk- Therefore, 
in the case when Afc = ak the eigenvalue problem (|3.1|1 is reduced to the equation (|3.6|l which 
is a first rank g-difference equation, i.e. 



(3-8) 

where Bq and ao are given by (|2.53p and (|2.54MTKK|l respectively. By applying the iteration 
method to (|3.8|) we find the solution 



iAg(l-g)g"xYao(9"~'=x) ' ^ ' 

where admissible choices of the real parameter ^k and function ^0 are presented in the table 
below (Figure Pi . A{x) is to be an arbitrary analytic function. 

Now, let us answer the question of when the solution t/'fc of (|3.9p belongs to the Hilbert 
space Tik- In order to do this we observe that 

I^.P eu) {x) = (1 ! ffjj^l^ Ok) {qx) . (3.10) 
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^0(2:) 


6 


7 > 












UQ — U 




-(7-l)A;-i 






6o = 6i = /i = 

7 / r\ 

62/0 




-(7- 1)A;- ilogq 


( b2-(l-g)A(0) 






7 = 




xA{x) 




7 < 




xA{x) 


-(7-l)A:-ilog,( 


^ 62-(l-q)A(0) 






62 = u 

/l / 

di an = ai 




-(7-l)fc-i 


1„„ Al-(l-qM{0)^ 




&2 = bi = /I = 
dian = ai 




-(7-l)A:-ilog,( 







Figure 1: Table of the forms of the function ^0 and the parameter 



The equation (|,S.10p can be written for 7 = in the form 

(iV'fcl^ Qk^ {x) = ^ ^|'^^^|2q^°^ (I^'^I^ Sk) (qx) , (3.11) 
and for 7 7^ in the form 

'\i^k\^ gk)ix)= (3.12) 

IV'fcP Qk) (qx) . 



j'-^b2{qx)^^ + bi{qxr + bo) ^, , ,2 



q'^ {b2 + llzf.Y'ZT"^ ) ('?"'^)'" + (1 - qrq-'-'Hq-'xp + bo 

We also observe that the function \^pk\'^ Qk does not depend on Aq{x). Using iteration method, 
after standard calculations we obtain the classes of solutions of (|3.1()p described in the following 
proposition. 

Proposition 1 For the solutions to the equation the following cases hold: 

1. For 7 = we have 

{\'^k? Qkj {x) = x"" , (3.13) 

where o"^ - g^fe+fei+^o) 
wnere q - (i_,)2q, • 

2. For 7 7^ we have following possibilities: 



(i) If bo / 0, 62 / and 62 + iiZ^^) ' .d,' ^ 0. then 

li^kf Qk) (x) = x'- \ \\' (°°),!_Vw^°^ • (3-14) 
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(ii) // 60 7^ 0, 62 7^ 0, /i 7^ and 62 + 



(l-g)^ (diao-ai) 
qdi 



0, then 



Qkj (x) = X 

(iii) Ifbo^O,b2^0,h = and 62 + '(\Z% ^"''7^7°'^ = 0, then 



V X2 



(iv) //60 = 0,h^0,b2^0, h^O and 62 + [IZ^^) ^''""°d7'^ ^ 0, i/jen 



(iV'fcl^ 0k^ (x) = x"" 



where q ^ = 



g2+7(fe-l)6i 



(v) If bo = 0, 61 / 0, 62 7^ 0, /i 7^ ant/ 62 + jiZ^^) ^'"'"^^r^^ = 0, then 



^2+7(fc-l)b-^ 



where q ^ 

(vi) //60 = = 0, 61 7^ 0, 62 7^ and 62 + jiZ^^) ^'''°rf7°'^ ^ 0' ^^^n 



{{'(pkf Qk) {x) = 



where q 



n>0; 



-9)^ (dion-oi) > 



X = a; 



where —q 



< 0. 



(vii) Ifbo = bi = 0,b2^0,h^O and 62 + [JZ^f) ^'"'7^;"'^ 7^ 0, then 



(a) 
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oo 



^h^re -q-'- = < 



(viii) Ifbo = bi = h = 0,b2^0 and 62 + \iZ% ^''^^"^^ / 0, then 



[m^Qk)ix)=x\ (3.23) 



where q 



gl-7+2fc7fc3 



1,, I (diap-ai) 
"2"^ (1-97) 9di 



/n a// the above cases xi,X2 are roots of the polynomial 

b2X^ + bix + bo = (3.24) 
and 2/1,2/2 are roots of the polynomial 

e'h2 + (1 - + (1 - ^)'^^-^/^- + ^0 = 0. (3.25) 

Proof: We easily obtain the subcases (i) — (iii) by iteration. The other cases are proved by 
calculation of the Laurent expression coefficient and application of Jacobi's identities 

k=—oo 

(see 0). 

□ 

The proposition given below classifies those function (|3.9|1 which are elements of Hilbert space 
'Hk- 

Proposition 2 The solution of equation i.V. 6]) belongs to the Hilbert space 7ik if and 

only if the parameters bo, bi, 62; a, h, di, uq, oi and j satisfy the following conditions: 

1. 1 = and 5;+^ < (1^. 



2. J > and one the following conditions is fulfilled: 
(i) 60/0, b2^0andb2+ llzf/'Z7'^ y^O; 



(ii) 60 / 0, 62 / 0, /i / and 62 + ;;:;^) ^^^';°,7"^^ = 0; 

(iii) 60 / 0, 62 / 0, /i = and 62 + {I'Zgf) ^''""gdT'^ = 0; 

(iv) 60 = 0, 61 / 0, 62 / 0, / 0, 62 + ^"'"X"^ / anrf A < 3^^^^; 

(v) 60 = 0, 61 / 0, 62 / 0, / 0, 62 + [iZ^f) ^'^X"^^ = anJ A < aZg!^ . 

(vi) 60 = = 0, 5i / 0, 62 / and 62 + {I'lgf) ^''^X"'^ ^ 0' 
(vii) in i/iis case the solutions never belong to the Hilbert space; 
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(viii) bo = h = h = 0,b,^0, b,+^i^^^ / and ^ < ^ (.^(^^^-^ - 1)) , 
The notation and classification given above are compatible with Proposition^ 
Proof: The function ipk belongs to the Hilbert space if 

j (l^fcl^ 0k) {x)dqx < +00 . (3.27) 

[a,b]q 

This is equivalent to 

oo 

5^(1 - q)q''y {li'lf ek) (q^y) < +00 (3.28) 

n=0 

for y = a,b. So, for the case (i) (i.e. 6o 7^ 0, 62 7^ and 62 + |^Zg7) ^'^^ grf^ 7^ 0) have 
from Proposition H that the jV-'fol^ Qk is given by H3.14jl . and we show that 



(l-g)2/^y ^l (^^7 ^- (^<+oo. (3.29) 



00 



From the identity 



where 



{q'''a-q')oo = jj;^. (3.30) 



{a;q^)^ = {l-a){l-q'a)... , (3.31) 
(a; g^)„ = (1 - a)(l - g^a) ... (1 - g^("-i)a) , (3.32) 



we obtain the conditions equivalent to l|3.29|l 



XI 



[^-q)y^T7^. — (°° , , x (3.33) 



Those conditions are fulfilled for 7 > 0. The proofs of the other cases are similar to the one 
above. 

□ 

Finally let us come back to the general situation and observe that (|3.4p . H3.5p and H3.6|l 
imply that the function 

rk{x) := A^ . . Al_^^,^pl_^{x) , n=l,...,k, (3.34) 

is an eigenvector of the operator with the eigenvalue 

Afc = dkdk-i ■ ■ ■ dk-n+iak-n (3.35) 

if i^k-n '■~ '^k-n is the eigenvector of Hk-n with eigenvalue ak-n- Moreover, one comes back 
to the eigensubspace C^};!_,„ acting on Cip^ by the annihilation operators Ak-n+i, ■■■ and 



13 



n 



(I • • • 




Figure 2: Presentation of action of the operators 



Afc. The above described procedures can be illustrated by a lattice of points in the {k,n) 
plane (Figure E]). 

The eigenfunctions of the operator Hfc given by (|3.9|1 and (|3.84p in the limit g — > 1 tend 

to 

V^^(x) = A^..AU+,x-('=-")(^-l)e-/o^/«W'^* for n = l,2,...k, (3.37) 
with the eigenvalues 

Xl = dk...di (^-aoik - n - 1) + ^{k - n) - b2-f^{k - n){k - n - 1)^ . (3.38) 

In the next sections we want to present some important examples, including the example 
of orthogonal polynomials of g-Hahn class which, in the limit (? — > 1, gives classical orthogonal 
polynomials. These examples will illustrate how the factorization method presented above 
works. 



4 g-Hahn orthogonal polynomials 

We obtain g-Hahn orthogonal polynomials when we require that the functions fk{x) = and 
dk = for /c E N U {0}. This is equivalent to 

7 = 1, (4.1) 

Bk{x) = Bo{x) = b2X^ + bix + bo , (4.2) 

Ao{x) = {[2]b2 - qao + q^ai) x + - (1 - q)h . (4.3) 

We see that the functions and are second and a first order polynomials respectively. 
From (|2.62|l we obtain that the function is also first order polynomial 

Akix) = q-^Aoiq-^x) + ^ Bq{x) = OfcX + b^ , (4.4) 
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where 



ak = {[2{k - l)]b2 + q-^ao - ai) , (4.5) 

h = -^-{l- q)q-''h . (4.6) 

Hence, the annihilation and creation operators are given by 

Ak = dg, (4.7) 

= -(52x2 + hx + bo)dqQ-^ - aux - h (4.8) 

and the Hamiltonian by 

Hk = - {h2X^ + bix + hQ)dqQ-^dq - (a^x + dq+ (4.9) 
+ 'Z''^'"'^ {-q-^a^[k - 1] + aiW - b2[k - l][k]) . 
The eigenvalue problem for the Hamiltonian (|4.9|1 is known as the g-Hahn equation [HI El 

{Bo{x)dqQ-'dg + Akix)dg) C = A^C • (4.10) 
The eigenvectors related to the eigenvalues 

AO = , (4.11) 

= akin] + b2[n][n - l]g-("-i) . (4.12) 

are given by 

4 = 1, (4.13) 

k 

t/;^ = A^..A^_„+il= H (^-(b2X^ + bix + bo)d,Q-^-a,x-b,)l, (4.14) 

i=k—n+l 

for A; G N U {0} and n = 1,2, . . . , k. The functions tp^ H4.14jl are polynomials. Each of the 
families {ipk}n=o ^ system of polynomials orthogonal with respect to the scalar product 
given by Jackson's integral 

ll^'^{x)ll^'^{x)Qk{x)dgX 5nm , (4.15) 

where the weight functions are obtained from (|2.68|1 

Bo{q '=+^x)...So(x) 

The classes of the weight functions go and the set of integration [a, b]g in (|4.15p are presented 
in HE]. 

In the limit g — > 1 this case gives us the classical orthogonal polynomials 

Bo{x)^ + Mx)4:) Pkix) = KPki^) ■ (4.17) 



dx^ dx 

The functions Bq and Af^ are second and first order polynomials given by 

Bk{x) = Bo{x) = b2x'^ + bix + bo , (4.18 
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Ak{x) = dkx + bk , (4.19) 

where 

ak = ~2{k - 1)62 + ai-ao, (4.20) 

bk = biih-k), (4.21) 

(in order to obtain this formulas we demand additionally that h = (^i^^^yA in l|4.6|l ). The 
eigenvectors tp'^ (orthogonal polynomials), in the limiting case, have the forms 

^l{x) = 1 , (4.22) 



^/.^(x) = (^o(^)^ + Ak{x)^ (^o(x)^ + Ak-i{x)^ . . . (^o(x)^ + Ak-n+i{x)^ 1 (4.23) 
and correspond to the eigenvalues 

XI = dku + b2n{n - 1) . (4.24) 

5 The case of constant weight functions 

We assume that all weight functions are constant Qkix) = const. We obtain two cases, which 
we consider below 

g-Deformation of the harmonic oscillator 

Additionally we demand that dk = i 60 = ^0 = 1 for the sake of transparency of the 
formulas. In this case we have: 

7 = 1, (5.1) 

Bkix) = 1 , (5.2) 

Ak{x) = , (5.3) 

fkix) = q-'foiq-'x) , (5.4) 

Qk = l , (5.5) 

where 



q'^jq ^ao~ai) h 1 1 1 a\ 

= V ^ (T^^ " (13^ • 

The annihilation and creation operators are given by 

Ak = dg + q-^Mq-'x) , (5.7) 
Al = -d,Q-'+q-'foiq-''x). (5.8) 

Solving the equation H3.8|l we find the basic state xp^ of the Hamiltonian given by 

Hfc = - (1 + (1 - q)q-'-'xfo{q-'-'x)) dgQ-'dg+ (5.9) 

+q~'' {Uq-'x) - q-'foiq-'-'x)) d,+ 

-q-'dgiUq-'-'x)) + q-''fi{q''x) + q-^^ (ao + {q^i - a^M) . 
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1. If ag 7^ qai, then 

V^O(x) = ^ \ — , (5.10) 

where xi and X2 are roots of the polynomial 

(1 - q)q^{q~^aQ - ai)x^ + (1 - qfhx +1 = (5.11) 

and G R\{0}. 

2. If ao = gai i /i / 0, then 

^fc(^) = ^ . (5.12) 
V(-(l - Q) hq ''x;q)^ 

It easy to see that the operator Q^^ acts as follows 

rO rO c? 1 c? , , , 

or equivalently 

ijlix) = g^Vo°(^) • (5.13) 
The functions are eigenvectors of the Hamiltonians H^. with the eigenvalues 

Xl = ak = q-^''{ao+iq^ai-ao)[k]) . (5.14) 

Similarly it is easy to obtain that the functions 

rkix) = Q-^r^ix) (5.15) 

are eigenvectors of with 

A^ = g-'M^o + (9V-ao)[fe]) , (5.16) 

in view the following commutation relations 

qAlQ-^ = Q-'Al_, , (5.17) 

AIQ = qQAl^, , (5.18) 

qAkQ-^ = Q-^Ak-i , (5.19) 

AkQ = qQ Ak+i . (5.20) 

Finally we present the action of the operators diagrammatically Figure 01 and state the fol- 
lowing 

Proposition 3 The functions 

C(x) = , ^ Q"-^A: . . . Al4ix) , (5.21) 

J{ao - gai)"ng!g"("-i)+'= 

for A; G N U {0} and n G N U {0}, where the function ipQ is given by 15. 1 or \5.1^) . are the 
eigenvectors of Hamiltonians \5.!J\) corresponding to the eigenvalues 

K = («o + {q^a^ - «o)[fc - n]) . (5.22) 
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n 




k 



Figure 3: Presentation of action of the operators 
In the hmit q ^ 1 this case gives us the harmonic oscillator 

Hfc = -^ + ^ X + + (ai - ao)fc . 

with eigenvectors 

V,n(^)=(^_^ + ^o_^^ye-^-^ for nGNU{0} 
corresponding to the eigenvalues 

= ao + (ao - ai)(n - k) . 



g-Deformation of the three— dimensional isotropic harmonic oscillator 

Additionally we demand that = and 6i = = 1- In this case we have 

7 = 2, 
Bkix) = 1 , 
Akix) = , 



(1 - q)x ' 



where 



2g^Q - ai) /l 



q^ (1 — q)x 
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The annihilation and creation operators have the form 

1 

Ak = d, + q-^'foiq-'x) - —\- , (5.32) 

(1 — q)x 

Al = -d,Q-' + q-''Mq-'x) - —\- , (5.33) 

(1 - q)x 

and the Hamiltonians are given by the formulas 

Hfc = - (q-" + (1 - q)q-"'-'xUq-''~'x)) d,Q~'d,+ (5.34) 



{ Ul-'x) - q-'Uq-'-'x)) d, - q~^^ { dMq'^-^x)) + + 



x^ 



+r''m-^x) + 2g-3'=M/o(^~^^) + q-^k (^a, + {q^ - ao)l^) , 

The basic states of the Hamiltonians H5.34p can be found as the solution (|3.8p . 
1. If ao / q^ai, then 

i^lix) = ' x^>^ , (5.35) 



where G M \ {0} and 



^k = -k + logq{l-q)Vh. (5.36) 
2. If ao = q'^ai, then 

V^^(x) = C°x^'= . (5.37) 
These are the eigenfunctions of the Hamiltonian corresponding to the eigenvalues 

Xl = a, = q-"' (ao + (/ai - ao)^) . (5.38) 

Finally we have the following lemma: 
Proposition 4 The functions 

V^^(x)=A^..A^,_„+lV^„= (5.39) 

/or n = 1, 2, . . . , fc, are i/ie eigenvectors of the Hamiltonian with the eigenvalues 

XI = q-'^-au-n = q-'^'"-''^ (ao + (q'a, - ^0)^^^^^^) • (5.40) 
In the limit ^ — > 1 this case gives us the three-dimensional isotropic harmonic oscillator 

dx'^ 16 
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ao-ai h 3ao + ai 
^('+2^ + ^^ 

with eigenvectors 



^Plix) = CfcOxt-'^e-^^' , (5.42) 



k / ~ \ 

d qq — ai hi i \ "o-ai ^: 



«W= n + + for „=l.^.t (6.43) 

i=k-n+l \ / 



corresponding to the eigenvalues 

XI = ao + {ai - ao){k - n) . (5.44) 

Appendix A. Derivation of formula (12.181) 

By the definition of the adjoint operator we have 



{Q*^kWk)k = {'4'k\QVk)k = J '^k{x)^k{qx)Qk{x)dqX = 

a 

oo oo 

J](l - (z)<7"feV^(?^V'fc(9"+'&)a(g"6) - - 9)'?"«V^^(^^fc(g"+'a)a('?"a) 

n=0 

oo 

^(1 - q)q'^bq-^J^{;^^^u{q'^b)Qk{q'^-H)- 



n=0 n=0 

oo 

m=n+l 



m=l 
oo 

- 5^(1 - g)(?'"ag-Vfc(9"^-ia)'/'fc(g'"«)a('Z""'a) ■ 

m=l 

In this sum the expression for m = 0, i.e. 

(1 - q) (bi^kiq~^b)ipk{b)Qk{b) - atljkiq~'^a)v^kia)Qkia)^ , (5.45) 

does not appear. The functions tpk{x) i ^k{x) are defined on the set {q^b : n e NU{0}}U{(7"a : 
n G NU {0}} and for the other points we shall put these functions equal to zero 

(g- V) (fe) = , 
(g- V) (a) = . 

From the equation H2.11|l we obtain for x / a and x ^b that 

h 

{Q*^kWk)k = j 'ipk{q~^x)ipk{x)Qk{x)^^^^^-^q'^dqX = ^g~-^ ^^-f^^ (Q~Vfc) 



k 

(5.46) 

Similarly we have 



{f'^kWk-l)k-l = j f{x)'>l)k{x)Lpk-l{x)Qk-l{x)dq 
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'4'k{x)f{x)ipk-i{x)rjk{x)gkix)dgX = {ipk\fVk^k~i)k , (5.47) 



where we use the equation (|2T 

Summarising we obtain the formula (j2.18p 

= (a, + fkT = Bk {-d,Q-^ + fk) -Ak{l + {1- q)xfk) , (5.48) 

where the operator is given by 

Appendix B. Derivation of formulas (12.221 - 127211 ) 

The operators of annihilation and creation given by H2.17( I2.18p can be rewritten in the form 

Afc = (9g + /fc = - ^ Q + ^k , (5.50) 
(1 — q)x 

Al = Bk {-dgQ-^ + fk) - ^fc (1 + (1 - q)xfk) = ^\ Q-^ + Vk^k , (5.51) 

(1 — q)x 

where the functions fk, % are defined by (|2.26p and H2.13|l . From the conditions H2.15MTT?)]) 
we have that 

rjk{qx)(pkiqx) = dkVk-i{x)fk-i{x) 

Bk{x)(pk{x) = dkBk-i{x)ipk-i{q~^x) (5.52) 
r]k{x)ipl{x) - dkr]k-i{x)ipl_-^{x) = dkUk-i - + g '^''^^-_^l^y]~J'''''"'^ . 

The first and second equations of (|5.52p are equivalent to 

Mix) _ ^^rjk^ 



iPk-i{x) Vkiqx) 

(Pkiqx) Bk-ijqx) 



iPk-i{x) Bkiqx) 
A simple calculation gives us 

Vk{x) = r]k-i{q~^x) = gk-i{x)rjk-i{q'^x) , (5.55) 

Bk-i{x) 

ipk{x) = dk ^^]^^^ fk-i{q'^x) = — ^^^v'fc_i(g"^x) , (5.56) 

Bk{x) gk-i{x) 

where the function gk{x) is given by H2.25p . Substituting H5.55| [5?56|l into the third relation 
in (|5.52p we obtain finally 

Vk-i{x)fl-i{x) - ^^^j^Vk-i{qx)ipl^i{qx) = (5.57) 

q^dkBk-i{qx) - gk-i{q'^x)Bk-i{q^x)\ gk-i{qx) 

"fcOfc-l — «A; H 



(1 — qYq^x"^ J d\ 
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